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Weak matrix elements: On the way to A/ = 1/2 rule and e'/e with 
staggered fermions. 

D. Pckurovsky and G. Kilcup^ 

'^Department of Physics, the Ohio State University, 
174 W. 18th Ave., Columbus OH 43210, USA 

We report progress in our study of hadronic weak matrix elements relevant for the AI ~ 1/2 rule and e' /e. The 
presented results are from our first runs on a quenched ensemble with f3 — 6.0, and a dynamical Nf — 2 ensemble 
with /3 = 5.7, using staggered gauge-invariant tadpole-improved fermionic operators. 



1. Introduction 

An important contribution of Lattice QCD to 
phenomenology is a first-principle calculation of 
non-perturbative matrix elements (ME's) of weak 
operators involving light hadrons. Knowledge 
of these ME's combined with experimental data 
translates directly to constraints on CKM matrix 
elements, and thus enables another test of the 
Minimal Standard Model. 

In this work we concentrate on weak decays of 
kaons into two pions. e' is the measure of di- 
rect CP violation in these decays. It is defined 
as e' = ilmA2e''^^^-^°'> where A2, 62, Aq 
and So are amplitudes and final interaction phases 
corresponding to isospin 2 and final pion state. 
New experiments will soon measure e' , hopefully 
resolving the currently muddy situation. In ad- 
dition to calculating e' , it is interesting whether 
Lattice QCD can adequately explain the domi- 
nance of AI — 1/2 transitions in kaon decays, i.e. 
the fact that lu = Re(AoM2) = 22.2. 

2. Method 

We work within the framework of an effective 
field theory obtained by integrating out W-boson 
and t-quark, using OPE and running down the 
effective Hamiltonian to scales of order of the lat- 
tice scale using RG equations @] . At these scales 
the effective Hamiltonian has the form 

r ^° r 1 



where r = -VtdVt*/VudV*^, Zi and yi are Wil- 
son coefficients (currently known at two-loop or- 
der), and Oi is the basis of 10 four- fermions oper- 
ators. In this effective theory e' can be expressed 
in terms of CKM matrix elements and {Tnr\Oi \K). 
Calculating the long-distance part of these ME's 
is the task for lattice theorists. 

It happens that the value of e' is dominated 
by competing contributions of the following two 
operators: 



Qa)v+A 



1 



9 



l{qi3(la)v+A 



(2) 



(3) 



In this talk we concentrate mostly on calculation 
of ME's involving Oq. 

In general, we follow the technique of calcu- 
lating ME's with staggered fermions developed 
in Ref. This method has proved to be suc- 
cessful in the past in calculation of the Bk pa- 
rameter, which enters the expression for indirect 
CP-violation parameter e. Shown in Fig. 1 are 
our latest results. During the last year, a new 
point at /3 = 6.4 has been added. The lattice 
spacings are determined by demanding asymp- 
totic scaling, with the overall scale taken from the 
continuum limit of the p mass. The final result is 
Bk{NDR, 2GeV) = 0.552 ± 0.007. 

We would like to achieve similar precision and 
finesse with much noisier ME's relevant for AI = 
1/2 rule and e'. We introduce a number of im- 
provements compared to the original work on 
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Figure 1. Quenched Bk using tadpole improved, 
gauge-invariant operators. 

these ME's [|. 

Due to technical complications, it is extremely 
difficult to extract four-point functions on the 
lattice. Instead, we calculate {tt\Hw\K) and 
{0\Hw\K) and (following Rcf. ||) use chiral per- 
turbation theory to relate them to {TnT\Hw\K). 

It is convenient to calculate B-ratios of ME's, 
which are defined as ratios of ME's to their val- 
ues obtained by vacuum saturation. We have 
to consider three types of 4-fermion contractions: 
'figure-eight', 'eye' and 'annihilation'. The latter 
two are notoriously noisy. See Ref. [Q for more 
details. 

3. Simulation 

Our simulation parameters are shown in Ta- 
ble |ll We use periodic boundary conditions, and 
replicate the lattice by a factor of 4 in the time 
direction. Only degenerate mesons {nis — nid — 
m„) are considered. We employ gauge-invariant, 
tadpole-improved staggered fermion operators. 

In Fig. 2 we show three contributions to the 
(bare) ME of Oq on a chiral plot. Although in- 
dividual contributions to Og might be diverging 
in chiral limit, their sum total seems to converge 
to a finite value with linear dependence on mq, 
in agreement with chiral symmetries of staggered 
fermions. 

4. Matching with continuum 

To quote continuum results we use the 'hor- 
izontal' matching procedure: using ffj^ as our 
renormalized coupling constant, we perform the 
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Figure 2. Chiral plot of ME involving Og operator 
for quenched ensemble, without perturbative cor- 
rections. Three contributions come from 'figure- 
eight', 'eye' and 'annihilation' contractions. 

matching to continuum at a scale q* with one- 
loop perturbation theory [^|j^]: 

-^C,j)Of +0(g4) + 0(a«) (4) 

and then run to a specified final scale 
/i (e.g. 2 GeV) using continuum two-loop equa- 
tions. We work in NDR variant of MS scheme. 

The finite coefficients Cy are uncomfortably 
big for some operators. In fact, for operators of 
type PP the perturbative correction is almost the 
size of the tree-level value. This makes one-loop 
perturbation theory unreliable. Directly related 
to this is a huge dependence on q* , which serves 
as an estimate of the size of second-order correc- 
tions. 

The situation can be made slightly better if we 
follow the continuum convention and quote re- 
sults in terms of a B parameter. While there are 
still potentially large perturbative corrections, we 
can minimize their effect by considering the B 
parameter for the operator Oe/Zp. This choice 
modifies both the matching and anomalous di- 
mension matrices, removing the largest correc- 
tions from the dominant PP matrix elements, at 
the expense of larger corrections in front of sub- 
dominant operators. As shown in figure 3, the 
remaining q* dependence is then modest. 

Of course, this sleight of hand doesn't solve 
the remaining problem, which is the apparent 
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Figure 3. ratio for quenched ensemble. The 
three groups of points correspond to bare opera- 
tors, and those renormalized in NDR scheme at 
2 GeV using two different values for q* . 

breakdown of perturbation theory for the pseu- 
doscalar renormalization Zp, or equivalcntly, the 
quark mass renormalization. Once this coefffcient 
is found (e.g. non-perturbatively) the value of 
will be readily available. The situation with other 
operators is the same. 

For the operators SS\ , 5^2 and PP\ the renor- 
malization coefficients are not yet known. For the 
central values we have therefore taken the coef- 
ficients appropriate for Landau gauge operators. 
Varying their values by 100% we find that the 
change in and i^s is insignificant. For other 
operators (namely and B'j) we will eventually 
need to determine the missing perturbative coef- 
ficients. 

5. Results and conclusions 

Our preliminary results for i?6(NDR, 2 GeV) 
at given values of /3 are 0.67 ± 0.04 ± 0.05 for 
quenched and 0.76 ±0.03 ±0.05 for dynamical en- 
sembles (first error is statistical, the second one is 
an estimate of higher-order perturbative correc- 
tions). For Bg/^(NDR, 2 GeV) we find 1.082 ± 
0.006±0.011 for quenched and 1.12 ± 0.02 ± 0.01 



for dynamical ensembles, in reasonable agreement 
with results using Landau gauge and smeared op- 
erators iQ. 

We do see a dominance of A/ = 1/2 tran- 
sition over A/ = 3/2 one. The ratio of am- 
plitudes varies sensitively with the quark mass, 
but at ma — .01, our preliminary results are: 
uj = Re^o/^2 = 11 ± 1 ± 3 for quenched and 
a; = 16±3±3 for dynamical ensemble. 

We have made a first step in a program to cal- 
culate e' and A/ =1/2 rule on the lattice. Apart 
from the common Lattice QCD problems of (par- 
tial) quenching, finite lattice spacing and size and 
degenerate quark masses, we have to face two 
more problems: failure of perturbation theory, 
and mostly unknown uncertainty in chiral per- 
turbation theory predictions. To solve the first 
problem, a non-perturbative determination of Zp 
is needed. We plan to repeat the calculations for 
a number of different /3's in order to take the con- 
tinuum limit. 

We thank the Columbia group for access to the 
dynamical gauge configurations, and the Ohio Su- 
percomputer Center for the necessary Cray-T3D 
time. 
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